REPORT No. 560 



A SIMPLIFIED APPLICATION OF THE METHOD OF OPERATORS TO THE 
CALCULATION OF DISTURBED MOTIONS OF AN AIRPLANE 

By Robert T. Jones 



SUMMARY 

A simplified treatment of the application of Heaviside's 
operational methods to problems of airplane dynamics is 
given. Certain graphical methods and logarithmic for- 
mulas that lessen the amount of computation involved are 
explained. 

The problem of representing a gust disturbance or 
control manipulation is taken up and it is pointed out 
that in certain cases arbitrary control manipulations may 
be dealt with as though they imposed specific constraints 
on the airplane, thus avoiding the necessity of any inte- 
gration whatever. 

The application of the calculations described in the text 
is illustrated by several examples chosen to show the use 
of the method* and the practicability of the graphical and 
logarithmic computations described. 

INTRODUCTION 

The theory of airplane dynamics in its present form 
is due mainly to the original researches of Lanchester 
and Bryan on the stability of airplanes.. Later inves- 
tigators, notably Bairstow and Wilson (reference 1), 
applied and extended the original conceptions of the 
theory. Bryant and Williams (reference 2) have re- 
cently shown how the operational mathematics of 
Heaviside may be used in applying the theory to prob- 
lems of the disturbed motions of airplanes. 

Although the calculation of disturbed motions of 
aircraft is important in problems of flight safety, little 
experience has been gained in the practical application 
of the theory owing to its mathematical complexity. 
The present paper gives the results of researches in the 
mathematical application of the theory. It has been 
found, as suggested by Bryant and Williams, that the 
Heaviside method affords the simplest and most direct 
solution of these problems. In order to bring out the 
advantages of this method, a treatment of its applica- 
tion is given and certain formulas and graphical con- 
structions are explained that make the calculations 
easier. 

In their usual form, problems of airplane dynamics 
depend for solution on the integration of simultaneous 
linear differential equations. Methods for the integra- 



tion of such equations are given by Wilson and Kouth 
(references 1 and 3) and in mathematical textbooks. 
The problems met in airplane dynamics are often more 
complex than the examples treated in textbooks and, 
when an attempt is made to apply the given methods 
to their solution, difficulties of computation arise. 

In view of the importance of investigating these 
problems and since their solutions involve lengthy cal- 
culations, it is desirable that as many mathematical 
simplifications as possible be employed. Heaviside's 
method gives such a simplification, the solution of ths 
differential equations being accomplished symbohcally 
by a single "expansion theorem." ' ^ 

THE DIFFERENTIAL EQUATIONS FOR THE DISTURBED 
MOTIONS • 

An airplane in uniform flight may be thought of as a 
free rigid-body in equilibrium. Deviations of the air- 
plane from this equifibrium condition may be caused 
by reactions due to control movement, gustiness in the 
air, or by some influence such as the stopping of an 
engine. The motions of the airplane following such a 
disturbance may be calculated if the momentary 
accelerations or forces are known. It is obvious that 
this computation may be performed by taking small 
intervals of the time and caculating the velocities and 
displacements generated by the known accelerations 
step by step, assuming the accelerations momentarily 
constant. 

The component linear and angular motions of the 
airplane in its deviations from equihbrium are given 
exact definition by constructing a set of axes rigidly 
fixed in the machine and considering its motions as 
being those of the axes themselves. The motions 
spoken of are then velocities and displacements of the 
airplane axes relative to the earth or the air. When 
the airplane is in steady flight, it maintains a certain 
equilibrium attitude with respect to the air and to 
the earth. Thus for climbing flight at a given engine 
speed a definite angle of attack and a definite angle of 
pitch must be preserved. Deviations from equihbrium 
in either sense will introduce reactions; hence motions 
of the airplane axes relative to both air and earth 
must be considered. 
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(A) "Lateral motions 

(B) Longitudinal motions — 



The aerodynamic reactions to the motions arise 
from changed relative air velocities over the different 
parts of the airplane. The calculation or measure- 
ment of these component aerodynamic reactions leads 
to quantities known as "resistance derivatives" or 
"stability derivatives," which are taken as constant 
factors of proportionality between the reactions and 
the velocities or displacements of the motions. For a 
more detailed exposition of the concept of stability 
derivatives, the reader is referred to standard text- 
books on aeronautics. 

On account of the bilateral symmetry of the airplane 
it is customary to divide the motions into two inde- 
pendent groups, the lateral and the longitudinal, each 
consisting of three degrees of freedom: 

EoUing. 
Yawing. 
Sideslipping. 
Htching. 

Vertical translation. 
Forward translation. 
Presumably the reactions to small increments of 
longitudinal speed or displacement do not sensibly 
influence the lateral motions and the two groups may 
be independently treated. In order to illustrate the 
calculation of the history of a motion due to a given 
disturbance, examples of lateral motions are chosen 
although the methods used are equally applicable to 
any set of degrees of freedom of the airplane. The 
quantities that arise in the consideration of the lateral 
motions are denned in the following table: 
Velocities and displacements of airplane axes: 
Z7 0 , equflibrium flight velocity along X axis. 
v, component of flight velocity along Y axis 

(sideslipping), 
p, component of angular velocity about X axis 
(rolling). 

r, component of angular velocity about Z axis 

(yawing). 
<p, angle of bank (relative to gravity). 
Forces and moments resolved along airplane axes: 
T, component of force along T axis. 
L, component of moment about X axis (rolling 
moment). 

N, component of moment about Z axis (yawing 
moment). 
Accelerations of airplane: 

Y 0 =Y/m (force per unit mass). 
io= L/rnkx 2 (moment per unit moment of inertia). 
N 0 =N/mk z 2 (moment per unit moment of inertia). 
Gust velocities resolved along airplane axes: 

Vo, component of gust velocity directed along Y 
axis. 

r 0 , component of angular velocity of gust about Z 
axis. 

p 0 , component of angular velocity of gust about X 
axis. 



Note. — The signs of the gust velocities are so 
chosen that a positive gust produces the same aero- 
dynamic reaction on the airplane as a positive velocity 
of the airplane in still air. The resolution of gust 
velocities along the moving axes is exact only t,o the 
first order of the small quantities involved. 

Airplane characteristics used as parameters: 

Y. 

Stability derivatives in terms of accelera- 
tions of airplane, thus : , 



. bYl 
bLI 
bNI 



A 



■ ^ 



Y T 

L P 
L. 
L r 
N, 
N„ 
N r 

"With the definition of the component motions that 
are to be considered, the stability derivatives will be 
of the form: 

bL bN bY 
3— > -y—> -5—1 etc. 
bp bv br 

where L, N, Y, respectively, are the rolling moment, 
the yawing moment, and the sidewise force, as they are 
customarily defined. 

It has been found convenient to transform all 
stability derivatives and disturbing effects into terms 
of accelerations of the airplane rather than retaining 
them as moments and forces. This transformation 
is accomplished by dividing out the appropriate 
moments of inertia and the mass of the machine. 

For example, ^jmk x 2 may be written simply as L p ; 

similarly jmk z 2 =N, and ^jp-j m = Y r - 

If the flight path is assumed to be horizontal (or 
nearly so) and the main forward velocity U 0 to be 
substantially constant, the equations of motion in 0 
lateral disturbance may be written: 

(In sideslipping) ^=9<p— rU 0 +vY,+rY r +Y 0 

(In rolling) ^=vL.+pL p +rL T +L 0 (1) 
dr 

(In yawing) ^- t =vN,+pN p +rN r +N 0 

In these equations the terms Y 0 , L 0 , and iV 0 represent 
known disturbing or controlling accelerations, assumed 
to be given as functions of the time t. In the first 
equation the terms g<p and —rU 0 are, respectively, the 
accelerations due to gravity and to the rotation of the 
moving axes. Since the axes chosen will ordinarily He 
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near the axes of the principal moments of inertia of 
the airplane, terms involving the products of inertia 
have been neglected. 

INTEGRATION OP EQUATIONS FOR VELOCITIES AND 
DISPLACEMENTS 

As previously mentioned, equations (1) may be inte- 
grated by taking small intervals of the time and calcu- 
lating the velocities, and finally the displacements, by 

dv 

assuming the accelerations ^> etc., to be momen- 
tarily constant. Although this method is sometimes 
useful, it naturally leads to extensive numerical work. 
The operational mathematics of Heaviside appear to 
offer the most promising means of performing these 
integrations. 

The first step in integrating the equations of motion 
by the operational method is to replace the symbol 
d 
dt 



by the so-called "differential operator" P, which is 



to be treated as though it were an ordinary algebraic 
quantity; the equations are then rearranged with the 
known disturbance effects on the right-hand, side: 

(D-Y c )v-g<p+(U 0 -Y r )r=Y 0 

-L,v+(P-L p )p-L r r=Io (2) 
-Nj>-N,p+(P-N r )r=N 0 

Since P<p=p, the. first equation may be operated on 
throughout by P, reducing all to the same variables 
(*» ft &• 

P(P-Y c )v-gp+D(U 0 -Y T )r=PY 0 

-L#+(P-L p )p-Ls=U (2a) 
-Nfi-N&+ (D-N T )r=N 0 

With the- equations in this form, they may be solved 
for v, p, or r by ordinary algebraic means; thus, 

DY 0 -g P{l\-Y T ) 

U (B-L,) -L, 



v= 



No 



(P-Nr) 



D(D—Y,) -g 
-L, (D-L p ) 
-N. -N P 



D(U 0 -Y r ) 
-L T 
(P-N T ) . 



(3) 



The expansion of the determinant of the numerator 
in terms of minors results in: 



(D-L p ) 
-N P 



(D—N r ) 



DY 0 + 



D(U 0 -Y r ) 

(P-Nr) 



In the calculation of any of the velocity components 
the same denominator appears; if this determinant is 
denoted by F(P), the forms of these components are: 



°—FlJD) Yo+ FiD) U+ TW) 0 



J 21 

V~FW) ^FW) ^^ 7 ^ ' 



MP)* 

FjD) 



(5) 



etc. 

Thus far the solution of the equations of motion has 
progressed simply on algebraic grounds, the required 
quantities (t>, p, etc.) having been found explicitly in 
terms of the symbol P. The symbol D was defined as 
the operation of derivation with respect to the time f, 
expressed by writing 

The terms of the solution f(D)/F(P) indicate that the 
formal operations are to be performed on whatever 
functions follow them as factors. Since they contain 
the symbol P in their denominators, it becomes neces- 
sary to define the operation indicated by 1/P or P' 1 . 
As P is an operation and not a number, its reciprocal is 
defined as the inverse of the operation of differentiation, 
rather than as the derivative itself divided into 1. 
The inverse of differentiation is integration: thus, 

D~ l =f . . . dt 

The operations indicated by the ratios of polynomials 
in P that occur in the terms of our solution then consist 
of a succession of differentiations [/(Z?)] and a succes- 
sion of integrations [F(D)]~ l . It is clear that the 
nature of the problems at hand requires that the result- 
ant of these operations be an integration, which is 
shown by the fact that the polynomial F(D) is invar- 
iably of higher degree in D than any of the polynomials 
f(P). 

THE EXPANSION EQUATION 

By treating the disturbances (such as Y 0 , N 0 , Lq) as 
discontinuous functions of the time, Heaviside obtained 
solutions of equations similar to the foregoing by a 
simple theorem. The substitution of Y 0 into Heavi- 
side's theorem results in 



KP) 

F(D) 



Y 0 =Y t 



/(0) 
W)' 



XF 7 ^)' 



(6) 



where the X's are the roots of the polynomial equation 
F(D)=0. This polynomial, F(P)=0, is used in the 
study of the stability of motion, being called the 
"stability equation." Its roots, Xi X 2 . . . X„, give an 



—9 
-N. 



-g D(U 0 -Yr) 
(D—L P ) -L T 



N 0 



(4) 
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indication of the natural tendencies of an airplane's 
motion and are used in the definition of stability. 

In order to apply the foregoing theorem to the 
integration of equations of airplane motion it is neces- 
sary to assume that the disturbance terms (Y 0 , N 0 , 
etc.) due to the control or gust in question are instantly 
applied at the assumed origin of the time scale (<=0) 
and remain constant thereafter. In the general case 
the disturbance terms in the equations of motion cannot 
be thus represented as remaining constant although in 
practical problems they may almost invariably be 
represented by means of functions of the form Y Q e nt . 
The interpretation of Heaviside's theorem (equation (6)) 
when this form of function is used is (see reference 2): 



/(X) 



(\-n)F'(X)' 



•1 



(7) 



"When dealing with variable disturbance terms, it is 
important to note that a discontinuity of the function 
representing the disturbance at t=Q is implied as in 
equation (6). 

By the substitution of (in) for n in equation (7), 
expressions that can be used when the disturbances 



■F(D)fbrJ)'X s tiy t 




Figure 1.— Map of polynomial. F (,D)=*I)i+<LD>+bD l +cI>+d near tero. 
jF(D)-»f (X) =0 when X>-*i±14. 

are represented by forms involving sin nt or cos nt are 
obtained 

-r{||> (cos „ (+iBb „ 0+ p^)|« «»]<7a) 

7 (ill) 

If i,,. ( =A+iB, then the expressions for the sine and 

-Z* \Vflt) 

cosine forms separately become 
(fD) 



F(D) 



Y Q sin nt 



= Y^A sin nt+B cos e><] 
f(D), 



(8) 



FTP) 



Y 0 cos nt 



--Y^A<M3nt-B8mnt+^(j^^ff^e^ (9) 



These latter forms are particularly useful because 
almost any arbitrary variation of gust or control may 
be expressed as a sum of sine or cosine terms. Thus if 

Y 0 =Ki sin nit+Kj sin «af + . . . + etc. 

TO 7o== Wo Kl 8111 nit+ iwj 1 8131 712 + 



. + etc. 
(10) 

Each of these terms may be evaluated by equation (8). 
SOLUTION OF OPERATIONAL EQUATIONS 

FINDING THE BOOTS OF THE EQUATION F(D)~0 

The expansion equations given for the forms 
f(P)/F(D) require the roots of the complementary 
equations F(JD).=0 for their _ solution. In cases of 
airplane motions this equation is normally of the fourth 
degree in D; hence it is not practicable to find the roots 
directly. Although a number of methods for approxi- 
mating the roots of such equations have been devised, 
the most direct way is to draw a curve of the function 
F(D) against D, locating the real roots as the points 
crossing the D axis. Usually in equations of this 
type near roots may be isolated by separating the 
equation into two parts. Thus, if 



F(D)=D i +aJD 3 +bD 2 +cD+d=0 



(11) 



there will usually be a large real root near D t =—aD 3 , 

or D=—a, and a small one near D=—~ This 

c 

division follows from the consideration that large roots 
are more dependent on the coefficients of the higher 
powers of D and small roots, on the lower powers. 

If the natural motion of the airplane contains oscilla- 
tory components, as it usually does, there will be pairs 
of conjugate complex roots of the polynomial F(D)=Q 
in addition to the real roots. The determination of 
these roots is naturally more difficult, although if real 
roots have been previously found they may be used 
to reduce the degree of the equation by synthetic divi- 
sion and the determination of further roots will be- 
come progressively easier. Complex roots of such an 
equation may be directly found by plotting a map of 
the polynomial F(D) for various values of D using the 
coordinates D=x+iy and finding the zero point, or 
root, by interpolation, as is shown in figure 1. If a 
very accurate value of the root is required it may be 
convenient to plot the region of F(D) near the origin 
to a magnified scale. Since the polynomial is what 
is known as an "analytic function" (reference 4), 



bF(D)_ 



i>F(D) 



(12) 



and the map in its smallest parts will consist of squares. 
In this way a more accurate interpolation may be 
made or a process analogous to Newton's method may 
be applied. 

It will be found most convenient to calculate the 
various values of F(D) by means of a vector diagram 
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as shown in figure 2. If trial values of D are expressed 
in the form B (cos 0+i sin 0) or Bc ie , vectors repre- 
senting each of the terms of the polynomial may be 
simply calculated. The problem is to make all terms 
of the polynomial balance each other and it is readily 
seen how this may be accomplished by varying 0 to 
change the relative inclinations of the vectors and by 
varying B to change their relative lengths. The 
advantage of this method is that it enables a close 
approximation of the value of a root with a minimum 
number of trials, the diagram making apparent how 
nearly all the vectors cancel each other. 

SOLUTION OF EXPANSION EQUATIONS 

The numerical operations indicated in the expansion 
equations (6) to (9) call for calculations with complex 
numbers (i. e., roots of F(D)=0). A great deal of the 
labor involved in these computations may be saved by 
the use of graphical and logarithmic methods. 

Thus, if it is desired to calculate values of the com- 
plex terms occurring in equation (6), the logarithmic 
formula 



/(X.) 



(13) 



log X#%y cM -M+log/(Xi)-log Xi-log F'ih) 

is used. For the purpose of calculating these loga- 
rithms, it is convenient to express the complex numbers 
(Xi, /(Xi), etc.) as vectors of radius B and angle 6, 
writing, for example, 

\ 1 =a+ib=B 1 (co3 Oi+i sin B^—B^h (14) 

by De Moivre's formula. 
A complex term of equation (6) may then be written, 

/(Xi) e x lt= 



~x7F(xD 





B^Bie"* 

B> . e Xif+f(9>-«i-«i> 

BiBz 



Then 



log^^^'=Xx«-|-log Bo+id 0 



(15) 



(16) 



and the resultant logarithm may be plotted as a straight 
line Xi*+ constant, which is then divided or extended to 
represent any division or extension of the time f over 
which the calculation is made. (See fig. 3.) The final 
vectors will represent the complex values of 

and it is seen that the ordinates of the points of the line 
Xif+constant are the angles of these final vectors while 
the abscissas are the logarithms of their radii. 



By a separation of the two components of the imagi- 
nary root \=a+ib, the logarithmic formula may be 
reduced to 



log 7 



/(o+fl) 



T «crt*» t = (a+ib) f +log Bo+id, 



'(a+ib)F'(a+ib) 

= (log Ro+af) +i(d 0 +bt) (17) 
The final formula, where X 1 =a+i6, then becomes 



/(Xi) 



or, by De Moivre's theorem, 



(18) 



/(X) 



XF'(X) 

The points 
spiral (fig. 3); 



e"=B 0 e at [coB (bt+6 0 )+i sin (bt+0 o )] (19) 



thus plotted will He on a logarithmic 
the deviation of this spiral from a circle 




/(#) 



Fiodee 2.— Graphical method of locating values of F(JD) near zero, where 
F(D) =ZM+o2)»+6X»+cX'+<i. Vi-n+iVi- R,c>'i-M<xa M-f sin «i) 

shows the influence of damping on the natural motion 
of the airplane. 

The summation indicated in equation (6) calls for 
the plotting of such a logarithmic spiral for each of the 
complex roots. Since these roots always occur in con- 
jugate pairs, the calculation may be carried out for one 
of such a pair and a spiral calculated for the second 
would be exactly conjugate to the first. Thus, it is only 
necessary to perform the foregoing calculations for one 
root of each pair, the summations indicated in the 
equations being carried out in effect by merely doubling 
the abscissas of the points of one of the conjugate 
spirals. If X L =a-H& and \ 2 =a—ib, this summation 
may be written: 

x» 

ZlxT^f xt =2£oe a 'cos Q>t+6 0 ) (20) 

The formulas for the integration of terms containing 
sin nt and cos nt may be put into a more convenient 
form for the graphical or logarithmic calculations, i. e., 

V 



j^Y 0 sin nt=Y 0 1 imaginary coordinate ^^/ n '+Zj ^+n^(\f 



n /(X) 



Y 0 cos nt- 



Y 0 real coordinate of 



/(fa), 
F(in) ( 



/(x) f 



(21) 



(22) 
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In these forms the graphical construction of the terms 

■^^e"" proceeds along the same lines as that of the 

terms involving complex roots X. Here the re- 
sulting diagrams -will he circles, divided into equal 
angles as ni may he divided. In case X is complex the 
plot of an e XI term will be a logarithmic spiral as before 
and it is important to remember that the summation 



For the velocities of an assumed gust, forms involv- 
ing e" are useful. Thus, if the gust is considered a 
"transient" one, disappearing rapidly from an arbi- 
trary initial value, the form (A) (fig. 4) 

t? 0 =flie-"' (24) 

may be used. Here Vi is the initial value and -n. is 
chosen to make the gust diminish in any required way. 




Figure 3.— Graphical construction ot £ ^l, ' x ' where \-a±Sb. Log /(X)-Iog X-]og F>(\)-\og Ro+lh. 



over each of a pair of conjugate roots is accomplished 
by doubling the abscissas of the spiral obtained for one. 

WAYS OF REPRESENTING GUSTS AND CONTROL 
MANIPULATIONS 

GUST DISTURBANCES 

If the disturbances to be considered are due to gusts, 
the terms Y 0 , L 0) etc. of equations (1) will be of the form 



Y 0 =v 0 Y t +r 0 Y, 



(23) 



where r 0 , p 0 , and r 0 are the component velocities of the 
gust, which may vary with the time. As given, the ref- 
erence system for specifying these gusts has been chosen 
so that a positive gust velocity may be considered as 
producing the same aerodynamic reaction on the air- 
plane as a positive airplane velocity in still air. All 
such gusts must be assumed to be moderate so that 
second-order effects 'may be neglected. (See refer- 
ence 1.) 



If the gust is to be made to start from an initial value 
of zero and to persist with the time, the form (B) 



(25) 



may be used. (See reference 1.) 

For the purpose of representing gusts that arise with 
any degree of sharpness from zero velocity to a given 
peak value and then diminish, the form (C) may be 
used: 

v 0 =Ke- nt (l-e- nt ) (26) 

The sharpness. of the rise of this gust is governed by 
-m and the decrease by -ra, since its curve approaches 
that of (l-e - "") near the origin and finally becomes 
asymptotic to e - "". 

In the- case of a rotating gust it is probably more 
logical to use -the transient forms that represent the 
gust as disappearing in time instead of being persistent. 

CONTROLLED MOTIONS 

When considering controlled motions, it is often just 
as reasonable to assume that the airplane is under a 
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kinematic constraint, or prescribed acceleration, im- 
posed by the control as to assume that the pilot uses 
the control in an arbitrary way. This assumption 
leads to the inversion of the integration problems here- 
tofore considered, because the motion of the airplane 
is itself predetermined and the forces and moments (or, 
more properly, accelerations) required to be supplied 
by the controls are calculated by differentiation. The 
ability of various control devices to produce a given 
maneuver of the airplane may thus be compared and 
the degree of coordination required of the-other-con- 
trols may be studied. 

The foregoing procedure is a particularly useful way 
of studying the lateral-control effectiveness in turns. 
Turn maneuvers, which usually begin and end in level 



naturally follows by differentiating this equation. If 
the turn is to be "perfect," that is, with no sideshpping, 
the rate of yawing throughout must bear a definite 
relation to the angle of bank, namely, 



-Biff erentiating the expressions for p and r gives the 
accelerations in rolling and yawing and hence the 




(C) v B - Ke- n \l-e- mi ) 
Fkjube 4.— Curves of different formulas for representing gusts. 

flight, may be described by means of a few sine or 
cosine terms. For example, the angle of bank <p may 
be given by 

<p=constant-r-^ii cos nt-\-Ai cos 2nt+etc. (27) 
(See fig. 5.) The rate of rolling at every instant 

130882—37 22 



or, simply, 



■ff sm <p 



r —jjV if v is under 30° 



(28) 
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Fioube 6.— Specifications for a turn maneuver In which the constraints are given by 
v-UiAi+A, cos n/+K4i cos 2nf; P-jf ; r— <?. 

moments, which will arise from two ' sources: the 
reactions due to natural stability and the reactions 
produced by the displaced controls. The reactions 
arising from the motions are found by combining the 
known stability derivatives with the angular velocities 
•p and r, obtained from the specification equations (27) 
and (28). The parts of the moments necessarily 
supplied by the controls are then obtained by deduct- 
ing these from the total moments. In the case of the 
aileron control, secondary moments in yaw result from 
tbe application of rolling moment, which modify the 
amount of rudder control displacement necessary. 

CONTROL AGAINST GUSTS OR ENGINE FAILURE 

In order to deal with attempted control of a given 
disturbance it is important to consider that there is 
invariably a lag in the pilot's reaction in countering 
the motion. In these cases it is possible to assume 
that the disturbance arises instantly, or nearly so 
(whether persistent or not), and that the pilot's dis- 
placement of the corrective control takes place accord- 
ing to the law 

5=50(1-6-"") (29) 

(see fig. 4(B)) where 5 0 is the assumed maximum con- 
trol deflection, which occurs more or less quickly as 
— n is made large or small. 
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EXAMPLES SHOWING APPLICATION OF OPERATIONAL 
METHODS TO PROBLEMS OF AIRPLANE MOTION. 

The following examples illustrate the application of 
the, various methods [- to specific problems of airplane 
motion. The airplane assumed in these calculations 
is a typical 2-passenger machine having the following 
characteristics: 

CHARACTERISTICS OF TYPICAL AIRPLANE 

Type: Monoplane; aspect ratio 6; rectangular, rounded 

tip, Clark T wing; dihedral angle, 1°. 
Dimensions: 

Gross weight 1,600 lb. 

Wing span 32 ft. 

Wing area 171 sq. ft. 

mk^ 1,216 slug-ft. a 

mk z * 1,700 slug-ft. 2 

Stability derivatives: 



Ob 
(fWsec) 


Cz. 




Zr 


Dot. 








160 
88.6 
63 


0.33 
L0 
'L8 


-S.44 

-3.23 
-2.46 


Lll 
1.88 
2.61 


-2.16 
—Lll 
-1.66 


-0.207 
-.301 
-.810 


-a 913 
-.663 
-.977 


6.82 
2.04 
L46 



» Flaps down. 

The calculated principal lateral-stability derivatives 
of this machine given with the other characteristics 
refer to motions of a set of axes fixed in the airplane 
but so inclined that the X axis points in the direction 
of the relative wind in straight flight at the lift coeffi- 
cient specified. The axes, nevertheless, move with the 
machine during the small oscillations considered and 
hence depart slightly from instantaneous reference axes 
fixed in the wind direction. 

ILLUSTRATION OF SOLUTION WITH CONSTANT DISTURBANCE 

TERM 

Example I, Rolling motion produced by deflecting 
ailerons at low speed: 

(a) Assume the machine to be in level steady flight 
at a speed of 88.5 feet per second ((7^=1.0) and that 
a rolling moment corresponding to <7 { =0.04, with an 
adverse yawing moment G n = — 0.01, is applied sud- 
denly at the time <=0. This condition corresponds 
approximately to a full deflection of ordinary ailerons 
at this speed. 

(b) The equations of the motion in the three degrees 
of lateral freedom may be set up without including the 
expressions for the lateral air force, since this force is 
Rrrm.11 and may be neglected in this case. The equa- 
tions are: 

dv 

T =gv-rTJ 0 



dr 



%- t =vN,+pN,+rN T +N 0 



(30) 



The terms L 0 and N 0 represent the accelerations due 
to the constant control moments suddenly applied at 
t=0. They are 

f- OaSb 



- iV ° Ink? 



(31) 



The substitution of D for d/dt, and the rearrange- 
ment of the equations result in 

m-gp+DUdr=0 

-La>+(P-L p ) P -L?=Lo (32) 
-Nfi-N„p+(D-N T )r=N 0 

Since the rolling motion is desired, the equations will 
be solved for p. The algebraic solution is: 



p= 



D 3 


0 


DU 0 


-L, 


Lo 


-L r 


-N. 


N 0 


(D-N r ) 


IP 


-9 


DU 0 


-L. 


(D-L,) 




-N, 


-N, 


{D-Nr) 



(33) 



which is then reduced to the form required for expan- 
sion in equation (6), 



(34) 



The calculation of the various polynomials in D results 
in: 



(35) 



f 1 (D)=L*-N T D i +U 0 NJ? 
f»{D)=L r D i -U 0 L t D 

F(D) =T>- (L P +N T )D 3 + (L„Nr-LrN p + U*N,)D* 
+ U 0 (L& P -I^.-%L,}D+g(LJfr-L r N,) 

At the assumed speed of 88.5 feet per second 
(0^=1.0), the constant rolling and yawing accelera- 
tions are 

1 ' 68 (36) 
N 0 = -0.301 

Using these numerical values, combined with those 
given for the stability derivatives, the polynomials in 
D become 

f 1 (D)=D 3 +0.66dD 2 +2.04D 

/ a (Z>)=1.88Z? J +l.llZ? (35a) 
J'(Z?)=Z> 4 +3.89Z? 3 +4.75Z) ! -1-10.33Z>— 1.13 

In order to perform the expansion of j^fcrj by Heavi- 

side's theorem it is necessary to determine the roots of 
the complementary equation F(D)=0. When the 
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polynomial F(D) is plotted as a function of a real vari- 
able (D)> two real roots of this equation are found: 

Xl= ~ 3 - 41 1 (37) 
X 3 = 0.104 J t3?) 

By the use of vector diagrams (see fig. 2) and the plot- 
ting of a map of the polynomial considered as a func- 
tion of a complex variable (D=x+iy), the following 
root was found by interpolation: 

X 3 =1.7S (cos 1.73-H sin 1.73) (38) 
An additional complex root that is the conjugate of X 3 
is known to exist and completes the four roots of the 
fourth-degree equation, 

X,=1.78 (cos 1.73— i sin 1.73) (39) 
The next step is to set up the integration equation 
and perform the indicated operations. Since the 
applied control moments Zo and N 0 are constants, form 
(6) will be used 



The logarithm of the result naturally occurs in the form 
x+iy. Plotting this point on the paper and construct- 
ing from it a line parallel to X3, we obtain the locus of 

10r X^F^ 

for various values of t (see fig. 3). The angles of 
the final points are given by the ordinates of these 
logarithms and the absolute lengths by the antiloga- 
rithms of the abscissas. The final points are found to 
lie on a logarithmic spiral whose radius decreases with 
the time (time measured as angle) showing the damping 
of this component of the motion. The summation 
over the two conjugate roots X 3 and Xj is accomplished 
without any further calculation by merely ddubling 
the abscissas of the points plotted above, as has been 
pointed out. The values thus obtained axe listed in 
the following table: 



\(D) MD) 



N 0 =Lo 



-0.484 results, 



=0.277 



(42) 



(43) 



The various terms to be substituted in this formula 
are found to be: 

F(0) U+ W) ° 

MX)2^+/ 2 (X)2V 0 =1.68X 3 +0.54X S +3.09X (41j 

X-F" (X) =4X 4 + 11.67X 3 + 9.49 X J + 10.33 X 
These terms are to be calculated for the four (real and 
complex) values of the roots. In the case of the real 
roots the calculation is made without resorting to 
graphical methods. For X!== — 3.41, the value 

/i(Xi)Z 0 +/ 2 (X 1 )iV 0 _ 

xTTO 

and for X a =0.104 

yi(X 2 )I 0 +/ 3 (X 2 )iV 0 

xlFtxl) 

It will be convenient to perform graphical calcula- 
tions to determine the other parts of the solution, 
corresponding to the complex terms. This result is 
accomplished by calculating the square, cube, and 
fourth power of the absolute length of X3 and by mul- 
tiplying each of these values by the proper coefficients 
in the polynomials /(D) and F(D ). By vector addition 
the value of the first polynomial was determined as 

/i(X 3 )io+/a(X 3 )iVo=3.99 (cos 5.23+i sin 5.23) (44) 
and the second 

X 3 i?"(X 3 )=40.6 (cos 5.60+i sin 5.60). (45) 
Since the quotient of these values is to be multiplied 
into e**' for a series of values of t it will be convenient 
to use the logarithm of this quotient, simply adding to 
it the various values of for which the calculation is 
to be made. This logarithm is 

X^'tXj) 




/i(0) 



Table of values obtained from graphical construction 



t 

(8econds) 




0 

.2 
.4 
.8 
.8 
L0 
LS 
2.0 


(ForXjaodXi) 


a is* 

.186 
.166 
.130 
.082 
-.030 
-.080 
-.110 



At the time tf=0, e x ' will be unity so that the initial 
condition of zero rate of rolling should be given by the 
sum of its coefficients. The summation 

-0.484+ 0.277 + 0.184=— 0.023 

shows how nearly this condition is attained. Figure 6 
shows the resultant rate of rolling and the components 
of the solution corresponding to each of the four roots, 
X„. In addition to the rolling curve obtained by the 
foregoing methods, other curves obtained by step-by- 
step integrations of the same equations of motion are 
given. In the calculation of these curves, steps of one- 
tenth and one-twentieth second were taken, which 
resulted in the differences shown. 

ILLUSTRATION OF SOLUTION WITH VARIABLE DISTURBANCE 
TEEMS 

Example II, Sideslipping during 2-control turn 
maneuver: 

(a) Assume the airplane to perform a specified bank- 
ing maneuver by application of a variable rolling mo- 
ment. If no yawing moments (from either rudder or 



log^ 



(log 3.99-log 40.6)+i(5.23-5.60) = -2.32-0.38 i 



(46) 
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ailerons) are applied, the natural stability of the air- 
plane will cause it to turn in a direction appropriate to 
the direction banked. Such a turn is called a "2-con- 
trol turn," inasmuch as only two (ailerons and elevator) 
of the three available controls are used. Since there 
will not be a very perfect coordination between the 
banking and yawing, some sideslip will result. It is of 
interest to know the approximate amount of this 
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Fiqdee 6.— Result of sample computation compared with step-by^etep Integrations; 
eiamplel. Bolllngmotlonfollowliigsuddenden^ctlonofanerons. Typical 1,600- 
ponnd airplane. Ci=L0; (7i»a04; C=— 0.01. 

sideslip during such a turn in studying the practica- 
bility of 2-control operation. 

(b) The first step in this problem is the determination 
of a suitable expression for the banking part of the 
.maneuver. It was considered that the pilot would 
naturally conform his use of the control to the desired 
motion of the airplane rather than move the control in 
a predetermined way and accept whatever motion of 
the machine followed. Hence it seems more logical 
to specify the banking motion itself rather than to try 
to predetermine a law of application of rolling moment. 

The airplane is thus assumed to be constrained in 
ban k in g by the aileron control so as to follow a well- 
executed bank maneuver and recovery. The usual 
procedure in making a turn is to bank the machine up 
to a definite angle, holding this angle steadily for a short 
time while in the steady part of the turn, and then to 
recover to level flight on the completion of the desired 
angle of turn. A curve representing such a relation of 
bank angle against time may be represented by a series 
of only two cosine terms with a constant defining the 



initial and terminal conditions of level flight, or zero 
bank angle. (See fig. 7.) For a fairly sharp turn with 
this small airplane the time required will be about 6 
seconds if the maximum angle of bank is 30°. The 
specification decided on is: 

Bank angle, p=0.327— 0.262 [cos t+% cos 2t] (47) 

which reaches a steady value of 30°, and gives level 
flight at t=0 and t=2ir seconds. The rate of rolling 
is the rate of change of this angle of bank; or 

#=^=0.262 sin f+0.131 sin 2t (48) 

A constraint of the machine in one of its degrees of 
freedom having thus been specified, it is only necessary 
to consider the equations for free motion in the remain- 
ing two degrees. As before, the lateral motion will be 
assumed to be independent of the longitudinal. There 
remain only the sideslipping and yawing motions to 
be considered. Their equations are: 



dv 
dr 

■fi=vN w +pN v +rN r 



(40) 



Although the equations contain the rate of rolling and 
the angle of bank, these are to be considered as known 
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FiouTUi 7.— Eesult of computation; example VL. Sideslip during a 2-control turn 

maneuver. 

from equations (20) and (30) and are, in fact, to be 
used as the disturbance terms. Calling 



(50) 



v g 
N 0 =pN, 

and rearranging the equations as in the other problems 



-UoN,^-+(D-N T )r=N a 



(51) 
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Solving algebraically for v/U 0 : 





T 0 1 




V 






u 0 


D 1 








-UoN, {V-N T ) 





(52) 



or 



where 



v_MD) Y ,MD) 



MD)=D-N r 

MD)=-i 

F(D)=D s -N T D+UoN, 



(52a) 



(53) 



If the airplane is to maintain its altitude while 
turning, the speed must be adjusted to give a higher 
lift than that at an equal lift coefficient in level flight. 
At an assumed lift coefficient of 1 the speed necessary 
to maintain altitude while turning at 30° bank is 
found to be 95 feet per second. Actually, if this 
speed is held throughout the specified maneuver, the 
longitudinal path will be accelerated somewhat; this 
condition will be neglected in the present problem. 
The necessary stability derivatives calculated for the 
new condition are: 

N r =- 0.712 
UoN,=2AO (54) 
N,= - 0.323 

The "disturbance effects" To and N 0 are (see equations 
(46), (47), and (49)) 



r o =0.111— 0.0888 cos t— 0.0222 cos 2t 
N 0 = -0.0846 sin t -0.0423 sin 2t 
and, finally, 

MD) 



(55) 



IT 0 - 111 " 



TOT 



cos t 



-0.0222^|^ cos 22- 0.0846^5} sin * (56) 
-0.0423^^- sin 21 

For the expansion of these terms in the integration 
equations (6), (8), and (9), it is necessary to know the 
roots of F(D) = 0. These are 



X=-0.356±1.51i 



(57) 



Since both these roots are complex, the operations 
indicated in the integration equations were performed 
graphically in the manner previously shown. 

The results of these calculations are shown in figure 7. 
The fact that the error in meeting the zero sideslip 
condition at the start of the maneuver was very small 
(even though the graphical construction of several 
terms was required) gives an indication of the accuracy 
of the calculation. 



Langlbt Memobial Aebonautical Laboratory, 
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Langley Field, Va., February 19, 19S6. 



APPENDIX 



EVALUATION OP ELEMENTARY OPERATORS 

A simple differential equation may be used to illus- 
trate briefly Heaviside's method of evaluating more 
elementary operational forms. Consider the case of an 
airplane executing pure rolling motion under the 
influence of a suddenly applied rolling moment of 
magnitude rnkxlq, which produces the impulsive accel- 
eration L 0 in roll. The equation of motion may be 
written: 

dt 



=pL p +L 0 



(58) 



in which both p and L 0 are supposed to have the value 
zero at the time t=Q. 

The solution of this equation as ordinarily found will 
consist of two parts, one of which is a solution of 



(59) 



the "complementary equation.'". In effect, Heaviside 
wrote both equations, (58) and (59), as one by intro- 
ducing a discontinuous function of t into (58). Thus, 
(substituting the usual D) 

Dp-L p ?=l(t)L<> (60) 

where the symbol l(t) is termed the "unit function," 
and is supposed to have the value zero until the time 
t=0 and to take the value 1 thereafter. The algebraic 
solution of (60) is then written 

p =D=r„ 1( - t)L ° (60a) 

and it is required to evaluate the form 

The procedure is to expand the fraction by the bi- 
nomial theorem in ascending powers of L p , thus, 

(D-L p )- 1 =D- 1 +D- 2 L P +D- 3 Lp*+ . . . +etc. (61) 

Since 



D- 1 l(t)=SHt)dt=l(t)t 



i 2 



D- 2 l(t)=ffl(t)dtdt=l(t)^ etc., 



(62) 



performing the indicated integrations results in 

{D-L p )-^{t)=l{t)(t . . . + etc.) (63) 

If this series is multiplied throughout by L P it becomes 
identically the series for e 1 * 4 except for the term 1, 
that is 

[^(D-L.y'+imt^Ht)^ (64) 



or 



(Z?-i )) )- 1 l(0 = ^(eV_i) (65) 
The final solution of the original equation (1) follows as 



p=l(t)^(e^-l) 



(66) 



Such forms as the left side of equation (60), involving 
the symbol D, are termed "operators." Equations 
(6) to (9) of the text are to be considered as evaluations 
of the more complex operators f(D)/F(P) along the 
above-indicated lines. The evaluation of a number 
of such forms is given in reference 5. 

Equation (6) of the text is a shorthand method of 
arriving at the foregoing solution. For the present 
problem this formula is: 



[Mr , YVoo-Ev, 



(67) 



and the various terms are: 



(68) 



f(D)=l 
F(D)=D-L„ 
/(0)=1 
F(0)=-L p 

/00=i 

F'(X) = 1 

The substitution of these terms in (67) results in 



as before. 



y=l(^(eV_l) 
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